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ABSTRACT 

We study the anisotropy of the cosmic microwave background (CMB) in 
cold and mixed dark matter (CDM and MDM) models, with non scale- invariant 
primordial power spectra (i.e. 1) and a late, sudden reionization of the 
intergalactic medium at redshift z r h- We test these models against recent 
detections of CMB anisotropy at large and intermediate angular scales. We 
find that current CMB anisotropy measurements cannot discriminate between 
CDM and MDM models. Our likelihood analysis indicates that models with 
blue power spectra (n ~ 1.2) and a reionization at z r h ~ 20 are most consistent 
with the anisotropy data considered here. Without reionization our analysis 
gives 1.0 < n < 1.26 (95% C.L.) for tt b = 0.05. 

Subject headings: cosmic microwave background - dark matter - galaxies: 
formation. 
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1. Introduction 

The COBE/DMR experiment ( Smoot et al. 1992 ) has revolutionized the field of 



structure formation, providing the first robust evidence for primary, large scale anisotropy 
of the cosmic microwave background (CMB). On these angular scales CMB anisotropies are 
expected because of potential fluctuations at last scattering ( |5achs fc Wolfe 19671) . These 



linear fluctuations are predicted to be constant in time, and then provide a unique tool for 
determining the initial conditions out of which large scale structure formed. Because of 
the low COBE/DMR angular resolution and of the correspondingly large cosmic variance 
( |Abbot fc Wise,1984|) , these initial conditions can not be yet determined with high precision 
( [Scaramella fc Vittorio 1993| , Bennet et al., 19941) . In spite of this, scale-invariant initial 



conditions, a robust prediction of inflation, are indeed consistent with the COBE/DMR four 
years data ( Pennet et al., 1996| ). It is worth remembering that the analysis of these data, 
based on the assumption of a flat universe, is completely insensitive to the chemistry of the 
universe. It seems difficult to reconcile the bulk of the observational data (CMB anisotropy, 
bulk flows, galaxy correlation function, etc.) with a model which is not dynamically 
dominated by cold dark matter (CDM). However, it does not seems probable that CDM 
can provide the closure density: if this were the case, the universe would probably be more 
inhomogeneous than observed on small scales. The excess power on these scales can be 
reduced, for example, by considering a low density (Qq ~ 0.2) CDM dominated model, 
which is flat because of a suitable cosmological constant, or a mixture of cold and hot dark 
matter, i.e. mixed dark matter (MDM) models. In this paper we will focus on the latter 
scenario. However, just because of the power reduction on small scales, in a MDM model 
structure formation tends to be a too recent process. This difficulty can be alleviated by 
considering initial conditions that are not exactly scale invariant. For example, it has been 
shown that the abundance of rich cluster in MDM models can be better reconciled with 
the X-ray luminosity function by assuming initial "blue" power spectra, i.e. P{k) = Ak n 
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with n > 1 ( Lucchin et al. 1996| ). Blue power spectra could overproduce small scale CMB 
anisotropy. However, a possible late reheating of the intergalactic medium could have 
controlled the level of anisotropy at or below the degree angular scale. If and when the 
universe underwent through a phase of early reionization of the intergalactic medium is an 
interesting and still open question, although from the Gunn and Peterson (|1965|) test we 
know that the universe must have been highly reionized at redshift z < 5, and the presence 
of heavy elements in the intracluster gas suggest that a considerable energy release occurred 
during the earliest stages of galaxy formation and evolution. The recent detections of 
cosmic microwave background (CMB) anisotropy at degree angular scales provide a new 
and powerful way for investigating this issue. In fact, the level of detected anisotropics 
at these scales can in principle discriminate among different assumptions for the thermal 
history of the universe. Also, using large and intermediate angular scales measurements 
increases the lever arm for determining the primordial spectral index. A definitive answer 
to these and other open questions will come when space missions (e.g. COBRAS/SAMBA 
and MAP) will provide a robust and definitive picture of the intermediate angular scale 
anisotropy. Meanwhile, the number of experiments reporting detections of anisotropy at 
degree angular scales has increased up to a couple of tens (see Table I below). Because of 
their sub-degree angular resolution, these experiments are quite effective in testing different 
reionization histories. So, the goal of this paper is to discuss if these anisotropy data, 
together with COBE/DMR, discriminate among different assumptions for the primordial 
spectral index n and for the thermal history of the universe. 



2. Theoretical calculations 

We consider a flat universe (f2 = 1) composed by baryons (0.03 < < 0.07), cold 
dark matter (Qcdm = 1 — — one family of massive neutrinos (Q u = 0.3), photons 
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and two families of massless neutrinos. For age considerations we fix H = 50km s~ 1 /Mpc, 
i.e. h = 0.5. The basic equations for describing the time evolution of density fluctuations in 
these different cosmic components have been already derived (Peebles 1981| , Peebles 1982a| , 
Peebles 1982b|) . In Fourier space, they are: 



i v + = $ (2) 

a 

/ i + ^f J i^f 1 V. T n e (v + 4^f 1 )-f / ^| (3) 

qJ a \qJ W l W v qJ) 

3 ,i / 3 - p! 2 - p 2 + 3p 2 p' 2 1 - /i /2 - 3p 2 (l - p' 2 ) \ I J T \ 
V ~ W-i \i_3//' 2 -/x 2 + 3/xV 2 3 - 3/i' 2 - 3p 2 (l - // 2 ) J \ Q J ^ 

ft, + 2-h = 8tcG (p B 5 B + PcdmScdm + 2p 7 5 7 + 2p v 5 v + 2p um lS. vm ) (4) 
\QttGcl 

tl 3 3-h= (p B V + Pumfvm + Pjf-y + Pvfv) (5) 

B= 2~ % ~^ 



v + H(t)v = a T n e ( / 7 - -v) (7) 

4dm = 2 (^) 

where $ = (1 — p 2 ) h — (1 — 3p 2 ) /i 33 . Eq.(l) describes the time evolution of the (Fourier 
transformed) fluctuation in the phase space distribution of massive neutrinos: bT = YQ, 
where Y = ye y (e y + 1)~ 2 , y = p/T v is the ratio between the neutrino momentum 



p = \/E 2 — m v and the neutrino temperature (we use natural units), proportional to 
the CMB temperature, T v = (4/H) 1 / 3 ^. The massive neutrino background density 
is then p um (t) = N um i:~ 2 T^(t)ip~ l where ip = / °° dyy 2 e (e y + 1)~ , N um is the number 
of neutrino families with a non vanishing rest mass (N vm = 1 in our case), e = E/T u , 
and each flavour state has two helicity states. We expand Q in Legendre polynomials, 
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Q = J2i=o9ePe(lJ'), an d transform Eq.(l) in a set of coupled differential equations for the 
multipoles g% ( see e.g. Ponometto et al. 1983|) . The massive neutrino density contrast 
reads <U(M) = ip' 1 J dyy 2 eYg while A um = ^ \ dyy 2 [y 2 + l/2(m„/T u ) 2 ]e~ 1 Yg . 
For e ~ y and e ~ m v , we have A um ~ 5 um and A um = 5 vm /2, respectively. Finally, 
Um = ip' 1 f dyy 3 Yg 1 . The integrals in y-space giving p um , A um , 8 vm and f um are performed 
numerically with a 16 points Gauss-Laguerre integration method. 

Eq.(2) describes massless neutrinos, and it is obtained by integrating in y- space 
Eq. (1) in the ultra-relativistic limit (e ~ y): l v = AQ. Again we expand l v in Legendre 
polynomials, l v = Y^e=o sePi(fi) and transform Eq. (2) in a set of 40 coupled differential 
equations. The density fluctuations of massless neutrinos is 5 U = sq, while /„ = si/3. 
We numerically follow the evolution of fluctuations in this component only when the 
perturbation proper wavelength is larger than one tenth of the horizon. Afterwards, free 
streaming rapidly damps fluctuations in this hot component. 

Eq.(3) describes fluctuations in the / and Q Stokes parameters of the CMB (see e.g. 
[Mclchiorri A. fc Vittorio, 1996| ) and differs from Eq. (2) because of the collisional term, 



proportional to the Thomson cross section, a T , and to the electronic number density, n e . 
The latter is evaluated assuming a standard recombination history QJones fc Wyse 1984|) , 
with a helium mass fraction Y# e = 0.23. We expand J 7 and Q 7 in Legendre polinomials: 
J 7 = J2i a iPi an d Q-y = J2i QiPi- The number of harmonics is automatically increased up to 
£ < 5000 for k < lMpc' 1 . 

The correlation function (acf ) of the temperature fluctuations can be written as 

1 oo 

C(a, B ) = — X^(2£ + l)C e P e (cas a) exp[-(£ + 1/2) 2 9 2 B ] (9) 

47F 1=2 

where 9b is the dispersion of a Gaussian approximating the angular response of the beam, 
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and 



A 2 A C OBE 



I 



dk k 2+n \a e (\k\)\ 2 , 



(10). 



32tt 2 {21 + l) 2 



We define the parameter A 2 = A/A<7 OB £ as the amplitude A of the power spectrum 
(considered as a free parameter) in units of Acobe, the amplitude needed to reproduce 
C(0°,4°.24) = (29fiK) 2 , as observed by COBE-DMR flBennet et al., 1996Q . 

Performing the integral in Eq.(10) with high accuracy at high £'s requires a very good 
sampling of the cr^'s in fc-space, and this is a heavy computational task. To avoid this 
problem we sample the interval —5 < logiok < with a step Alog^k = 0.01, we evaluate 
the integral in Eq.(10), and we use a smoothing algorithm to suppress the high frequency, 
sample noise. The Cgs obtained in this way differ from those obtained with a much denser 
/c-space sampling (Alogiok = 0.001) by only a fraction of percent up to £ < 2000. 

The time evolution of the baryon and CDM density contrasts and of the baryon 
peculiar velocity are described by Eq.(6), (8) and (7) respectively. The system is closed by 
Eq.(4) and (5) describing the field equations for the trace and the 3 — 3 component of the 
metric perturbation tensor. 

We numerically integrate the previous equations from redshift z = 10 7 up to the 
present. In the following we will also assume that the universe reionized instantaneously 
at redshift z r h « 1000, and remained completely reionized up to the present. Because of 
the low baryonic abundance, we need reionization at high redshifts in order to substantially 
suppress the anisotropy at degree scales. For Q b = 0.05 and z rh < 70, there is a probability 
lower than 50 % for a photon to scatter against a free electron at z < z rh . 



- 8- 



3. Data analysis 

We have selected a set of 20 different anisotropy detections obtained by different 
experiments, or by the same experiment with different window functions and/or at 
different frequencies. For each detection, labeled by the index j, we report in Table 1 
the detected mean square anisotropy, A^ xp \ and the corresponding \-a error, Ylf xp \ 
In our error estimates, the calibration error was added in quadrature to the statistical 
error. When not explicitly given, we estimated the mean square anisotropy as follows: 
^(exp) _ y ar [{/¥p.y^ _ a i where {ATj}j are the published anisotropy data of the j-th 
experiment, l^ar^ATj},,] is the variance of the data points, and <j| is the mean square value 
of the instrumental noise. 

Theoretically, the mean (over the ensemble) squared anisotropy is given by a weighted 
sum of the C/s: 



[ A i] = ^ E(2* + l)C e Wu = £(2£ + l)W ej (13) 
where the windows function Wgj contains all the experimental details (chop, modulation, 



beam ,etc), and Cp eff is the mean value of the C/s over the window function. The effective 
multipole number £ e ff,j is defined as follows: 

and is listed in Table I for a scale invariant model without reionization. Although in 
principle model dependent, the values of £ e ff,j are quite stable because of the narrowness in 
£ space of the window functions. 

Using numerical simulations, which take into account scan strategy and experimental 
noise, we verify that the expected distribution for Af xp ^ is well approximated by a gaussian, 
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with mean and a cosmic/sampling variance 



zr i) = y^U2l + l)Wlf} (14) 

Here fj represents the fraction of the sky sampled by each experiment and it is also listed 
in Table I. 

Given the Gaussian distribution of A^ xp \ we compute the likelihood of the 20 (assumed 
independent) CMB anisotropy detections as follows: 

1 i [A (ex ' p) - A 2 A {th) ] 2 

c(An ' zrh) ^ jmw^wr { ~ 2 ^Af¥ + ^ } (15) 

As already stated, this is a function of three parameters: the amplitude A, the spectral 
index n and the reionization redshift z r h- 

For each pair n — z r h we select the value A m ax which maximizes the Likelihood 
(isolevels of A max are shown in Fig.l for Q b = 0.03,0.05,0.07, respectively). The results 
of our theoretical calculation are shown in Fig. 2a, where we plot the C/s for few models, 
each of them normalized with its own A ma x- Our results show that the difference between 
a pure CDM and a MDM (Q u = 0.3) models, normalized with the same A max is very tiny, 
< 2% up to i < 300 and < 8% up to i < 800 ( |de Gasperis et al., 1995| , |Gates, Dodeiioil 



and Stebbins 19951 , |Ma fc Bertschinger 1995Q . So we will make hereafter the assumption 



that the anisotropy pattern does not depend on Q u , at least for Q u < 0.3. This generalizes 
our assumption that the CMB anisotropy depends only upon three independent variables: 
A, n and z r h- For the models shown in Fig. la we also plot the band-power estimate for the 
j-th experiment Ci eff = Air / (2£ + l)Wij and in Fig. lb the corresponding window 
functions, Wi t j. 

The 2-D, conditional distribution L(n, z r \A max ) = C(A max , n, z r h) has a quite distinctive 
peak at n = 1.24 and z r h = 20. A P-confidence level contour in the n — z r h is obtained 
cutting the L distribution with the isolevel Lp, and by requiring that the volume below the 
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surface inside L P is a fraction P of the total volume. In Fig.3 we plot the L 68 and the L 95 
contours, again for Q b = 0.03, 0.05 and 0.07, respectively (we actually sampled a region of 
the n — z r h plane much larger than shown in Figure). 

It is clear that the considered data set identifies a preferred region of the n — z rh space. 
The shape of the confidence contours can be understood by noting that increasing n would 
overproduce anisotropies: a corresponding increase in z r h is thus required to damp the 
fluctuations to a level compatible with observations. On the other hand, if f2 b is increased, 
the optical depth increases, thus requiring a lower reionization redshift to produce the 
detected level of degree-scale anisotropy. This effect dominates over the smaller increase of 
primary CMB anisotropies. 

The simple analysis carried out here does not take into account details on the scan 
pattern and/or the beam profile. Moreover, we know that the published error bars could 
not account properly for correlated errors or other more subtle effects. In order to test 
the significance of our analysis and its robustness against the published estimates for 
the experimental errors, we take a drastic point of view: we consider the best fit model 
(n = 1.24, z rh = 20) normalized to COBE/DMR (i.e. we fix A = 1). For each degree-scale 
data point we compute the difference between the measured value and the value expected 
from the theory. We then associate to each data point a 1-a error bar equal to the root 
mean square of these 19 differences (instead of the published error bars). The corresponding 
results of the likelihood analysis are in very good agreement with what is found using 
the appropriate errors: n > 1 is slightly better, especially if the reionization redshift is 
increased. So we do believe that our results are little affected by possible correlations in the 
errors and/or systematics. 
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4. Conclusions 

Our main conclusions, derived from Fig.3, are as follows. 

1. The conditional Likelihood shows a maximum for n = 1.24 and z r h = 20, for 
0.03 < Qb < 0.07, i.e. for baryonic abundances consistent with primordial nucleosynthesis. 
We have checked the stability of this result by applying a jack-knife analysis to the 
considered data set, and we have seen that the results are reasonably stable: the best model 
has always n = 1.24(±0.02) and z r h = 20(±10), unless COBE data are excluded (in that 
case n = 1.40 and z r u = 20). 

2. The 95% confidence contours in the n — z r h plane include a wide range of parameters 
combinations. This means that the presently available data set is not sensitive enough to 
produce "precise" determinations for n and z r h\ systematic and statistical errors in the 
different experiments are still significant. 

3. If we exclude an early reionization of the intergalactic medium (z r h = 0) we get 
the following 95% confidence level estimates for the spectral index: 1.02 < n < 1.28 
(Q b = 0.03); 1.00 < n < 1.26 (fi 6 = 0.05); 0.96 < n < 1.24 (ft b = 0.07). This has to be 
compared to the results from COBE-DMR alone: n = (1.3 ± 0.3), at the 68% confidence 
level ( [Bennet et al., 199~6| ). So, in spite of their still low signal to noise ratio, the degree 



scale experiments already allow to better constrain the spectral index, although still at 
the 10% level. Note that the "standard", flat model with no reionization is close but not 
always inside the 95% confidence contour. However, this result is determined mainly by 
the Saskatoon experiment. In fact, we tested the stability of our results by repeating 
the analysis several times, excluding one experiment at the time. If we do not consider 
reionization, the 95% C.L. interval 1.00 < n < 1.26 does not change more than a few % 
excluding either Argo, CAT, MAX, MSAM, South Pole or Tenerife. The lower limit drops 
to 0.86 < n if the Saskatoon data are excluded, while the upper limit raises to n < 1.76 if 
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the COBE data are excluded. This test also shows that neglecting the correlation due to 
overlapping sky coverage (e.g. Tenerife and COBE, and/or MSAM and Saskatoon) does 
not change significantly the results of our analysis. 

4. For scale invariant models with no reionization, the height of the first "Doppler" 
peak occurs at £ ~ 220 and is a factor of ~ 5.6 higher than the Sachs- Wolfe plateau at low 
£'s. For n = 1.2 the peak amplitude is roughly a factor 1.5 higher than in the scale-invariant 
case. A complete reionization from z r h — 20 up to the present suppress the peak by roughly 
20%. Altogether, a model with n = 1.2 and z r h = 20 has a Doppler peak a factor of 2 higher 
than in the standard scale invariant case without reheating. So our analysis confirms that a 
Doppler peak in the C% spectrum centered at ~ 200 is perfectly consistent with the data. 

All the models we have worked with have, as stated above, h = 0.5. In flat models 
with vanishing cosmological constant we can expect 0.4 < h < 0.65 from the estimated 
age of globular clusters ( [Kolb E.W. fc Turner JVLS1, ). Varying h modifies the amplitude 



and position of the acoustic peaks in the radiation power spectrum. However, we have 
verified that, in the region of the spectrum probed by current anisotropy experiments, 20% 
variations of the Hubble constant yeld modifications in the spectrum very similar to those 
obtained varying Qf, between 0.03 and 0.07. Also, we did not consider tensor modes in 
our analysis, as they are expected to be of negligible amplitude for n > 1 flSteinhardt P.J 
1995|) . We will address this issue in a forthcoming paper. 
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Fig. 1. — Isolevels of the normalization amplitude A max , which maximizes the combined 
likelihood £(A, n, z rh ). 

Fig. 2. — Power spectra of CMB anisotropies for different combinations of spectral index 
and reionization redshift (panel a). The data points are derived from the experiments listed 
in Table 1. In panel b we plot the corresponding filter functions. 

Fig. 3. — Confidence level (68 and 95%) regions for the spectral index n and the reionization 
redshift z r h- The black square marks the position of the model featuring the maximum 
likelihood. 



TABLE 1 

CMB Anisotropy Detections Used in the Analysis 



experiment 


ref. 


AT 2 (fiK 2 ) 


68%+, (iiK 2 ) 


68%-, (iiK 2 ) 


sky coverage 


*eff 


COBE 


1 


841 


58 


57 


0.6580 


5.8 


Tenerife 


2 


1770 


840 


670 


0.0124 


20.1 


South Pole Q 


3 


480 


470 


160 


0.0005 


49.4 


ooutn r oie iv 


A 
4 


on/in 


OQQfl 

zooU 


< yu 


U.UUUO 


00. / 


Python 


5 


1940 


1890 


490 


0.0006 


129.0 


ARGO Hercules 


6 


360 


170 


140 


0.0024 


118.9 


ARGO Aries 


7 


580 


150 


130 


0.0055 


118.9 


MAX HR 


8 


2430 


1850 


1020 


0.0002 


162.0 


MAX PH 


9 


5960 


5080 


2190 


0.0002 


162.0 


MAX GUM 


10 


6580 


4450 


2320 


0.0002 


162.0 


MAX ID 


11 


4960 


5690 


2330 


0.0002 


162.0 


MAX SH 


12 


5740 


6280 


2900 


0.0002 


162.0 


MSAM 3F 


13 


4680 


4200 


2450 


0.0007 


181.3 


Saskatoon 


14 


1990 


950 


630 


0.0037 


99.9 


Saskatoon 


15 


4490 


1690 


1360 


0.0037 


175.4 


Saskatoon 


16 


6930 


2770 


2140 


0.0037 


235.2 


Saskatoon 


17 


6980 


3030 


2310 


0.0037 


283.2 


Saskatoon 


18 


4730 


3380 


3190 


0.0037 


347.8 


CAT 1 


19 


1180 


720 


520 


0.0001 


414.9 


CAT 2 


20 


760 


760 


360 


0.0001 


579.7 
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1996; (13) Cheng et al. 1994; (14, 15, 16, 17, 18) Netterfield et al. 1996; (19, 20) Scott et al. 1996 
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